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Abstract: The instability of vacuum energy in de Sitter space as discussed recently 
by Polyakov is argued to be a generic feature when external gravitational fields are 
present. It is related to the existence of forbidden (by momentum conservation) 
decays derived in some detail by Bros, Epstein and Moschella. Some calculations are 
discussed in a conformally invariant setting, and a fiat space model is presented. 
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1. Introduction 

The understanding of the cosmological constant and its relationship with the vacuum 
energy is one of the most important problems in theoretical physics as of today. A 
recent contribution to it is Polyakov's [1] claim that vacuum energy is unstable. This 
instability is supposed to be different from previous proposals in the sense that this 
is due to particle interactions in a de Sitter background, and consequently is only 
indirectly of gravitational origin. 

To be specific, what is asserted is that interactions draw a generic instability as 
an imaginary part of the free energy of the quantum fields as computed in a de Sitter 
background. 

As a matter of fact Polyakov's is a three-loop effect and, as such, lies outside of 
our previous analysis [2], which was restricted to one-loop order. Our purpose in the 
present paper is to briefly comment on some physical aspects of this effect. 

To begin with, it is worth pointing out the following. In an interesting series of 
papers, Bros, Epstein and Moschella [3] following early work ^ in [4] and [5], have 
shown that one particle decays in (f)^ or (j)^ theories are not forbidden kinematically 

^In the book by BirrcU and Davies [6] some earlier references can be found. 
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in de Sitter space. Representing by h{x) the scalar field, such decays imply a nonva- 
nishing width 

r {h — > hh) or else F {h hhh) 

This is in sharp contrast with the situation in flat space, where momentum 
conservation forbids them. The reason for that is the lack of translational invariance 
(there is no abelian translation subgroup of the de Sitter group; pseudotranslations 
do not commute), so that two-point functions are not necessarily functions of the 
difference between spacetime coordinates of the two points, which is the root of global 
momentum conservation in any physical process. In fact this effect is common to 
any quantum field theory in a nontrivial gravitational background. 

Once momentum conservation is not working, nothing forbids the vacuum de- 
cay in to physical particles, which essentially related to the effect pointed out by 
Polyakov. Assuming, as we do, crossing symmetry, the preceding channels are re- 
lated to the vacuum decay in the tree approximation 

r (0 ^ hhh) or else T (0 ^ hhhh) 

Once there is a nonvanishing amplitude for this sort of decay into several parti- 
cles, it seems plain that the inverse reaction is much less likely, so that there is an 
enhanced production until the particle density n is so high that 

^ r 

at which point detailed balance should establish itself and the particle production 
growing stops. 

In a curved space the concept of particle depends on the observer through the 
definition of positive frequencies associated to a complete system of solutions of 
the classical field equations. For example, if fl{x) and ff{x) are two such systems 
appropriate in the far past and in the far future of the process of interest, the free 
quantum field can be expanded as 

ct>{x) = J2 («n/^(a:) + aifrix)) = J2 {Kflix) + hifl*{x)) 

n n 

This defines the operators a„ and 6„. This is also true, in particular, for the asymp- 
totic fields, and 0out- Now four different vacua can be defined as usual in an 
S-matrix approach 

<m = aTKj = 

The two sets of modes are related by a Bogohubov transformation, so that a 
given vacuum contains particles as defined with respect to a different vacuum (that 
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is, using a different definition of positive frequency). Wlien interactions are taken 
into account even with the same definition of positive frequency the vacuum is not 
stable 



In practice botli effects (tfiat is, particle creation due to the external gravita- 
tional field at zero coupling as reported, for example, in [7], and the effects of the 
interaction) compete, and in order to separate them one has to study specific chan- 
nels as well as their dependence on the coupling. 

The physical matrix elements to calculate in a a — >■ /3 are referred to the ap- 
propriate notion of particles out(/3^|«^)in, i-e. choosing the modes corresponding to 
the (asymptotic) past and future zones where the field behaves as a free field. Even 
when S-matrix elements are not defined sensu stricto (such as in de Sitter space) 
transition amplitudes for finite time intervals can still be computed using Feynman's 
rules. Our point of view is similar to the one in [8] in that we assume that enough 
of the analytical scheme of fiat space quantum field theory survives to justify the 
formal use of the interaction representation and related path integral techniques. 

In order to separate the true interaction effects from the creation of particles due 
to the gravitational field, we need to focus in the out(/3|«)in matrix elements, referred 
to a common notion of vacuum, which are related by a Bogoliubov transformation 
to the previous ones. 

For this matrix elements we can define the S'-matrix as 5*^10)^ = |/3)out, which 
has the familiar interaction representation form 



In any ^-matrix perturbative framework, unitarity precisely relates the imagi- 
nary part of the vacuum diagrams to creation and absorption of physical particles 
from the vacuum, 



Unitarity means that for any couple of states \a) and and any closure relation 



l(o^;ft|o^/)l^i 




s = i + ir 




the following is true 



{a\h) = {a\SS^h) = (a| (1 + iT) (l - iV) \h) 



so that 



{a\i {T-V) \h) = -(alTT^Ife) = -^(a|r|n)(n|rt|6) 



n 



3 



In a A0'^ theory, to second order we have for the vacuum-to-vacuum amphtude: 

{o\i (r(') - r(2)t) |o) = _ j](o|r(^V) • (n|r«^|o) 

n 

where we do not specify the appropriate asymptotic hmit (i.e. in or out), since the 
obtained result is identical 

2im(o|r(')|o) = ^l(o|r«|n)r (1-1) 

n 

Up to second order, the T matrix for a X4>'^ theory is 
S = Texp J dx : (l){xf :| 

y 1 j dxdyT{:4^{xY:x:(t>{yY ■■) + ■■■ (1-2) 



4! 



so the previous relation gives^ 



Re j dxdyG{x,yf^ j dxdyW{x,yf (1.3) 



where G{x,y) = (O|T0(x)0(y)|O) is the "Feynman" propagator and W{x,y) 
{0\(f){x)(f){y)\0) the Wightman function. 



2. Vacuum decay 

In Poincare coordinates the metric of de Sitter space reads 

ds^^ (^^y {du^ - d^') 

A conformally coupled scalar field is massless and the value of the curvature 
coupling is ^ = l^^l' where n is the dimension of spacetime. In de Sitter space, 
where the curvature is constant, R — "^"^^^ , this is equivalent to a minimally coupled 
(C = 0) scalar field with mass = = 

This mass is in the complementary scries of the dc Sitter group SO{n, 1), with 
a,iji = ^ parameter (m^/^ = ji^ + ^^^l^ ) ■, and the functional form of their two-point 
functions (without entering into the ie prescriptions for the time being) is particularly 
simple 

r(f-l) p /n n _ n. l + ^ T (| - l) _ _„ 
Z-2(47r)t V2'2 '2' 2 2(27r)tZ"-2 ^> 

^We ignore the tadpole contributions that disappear by considering normal order in the inter- 
action term. 
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Let us concentrate in this simplest example for the time being. In the coordinates 
we are using 

1 ^ _ -{u-u'f + {x-x'f 
'~ 2uu' 

Exploiting the conformal invariance of the setup ^, we can expand the free field 
expansion in term of the ra- dimensional modes 

n/2-l 

/k(M,x) = ^e-*fc«gikx 

^ ' (27r)^K2-i^ 

which are equivalent to choosing the Euclidean or Bunch-Davies vacuum. 



It is a fact that the closure relation for the complete set of solutions gives a 
particular solution of the homogeneous equation (id est, without the delta function 
source) namely 



/ 



d"-^k/k(n,x)/*(u',x') = ^^^^^ / ^""^k ^,uiu-u')^H.-.') 

r(f-l) / 2uu' 



2(27r)2/"-2 \{t{u - u' - ie)y + 
where the integrals can be done with 

^ n — 3 

dy{l - yr-^e-y = (0 r - l) J^(a 



(2.2) 



^ 0F(«2 + ^2)t-i 

and the last one needs an small real part for a. This closure relation appears in 
the derivation of (1.3), and this particular homogeneous solution is precisely the 
Wightman function. It is plain in (2.2) that the ie prescription depends on the time 
ordering of the arguments. The unitarity relations rely on this subtle difference with 
the Feynman propagator. 

In the X(j)^ theory, the bubble diagram G{x, yY can be computed for the con- 
formally coupled case in 4 dimensions. The propagator is — where z 

z — \ — ie 

is the geodesic distance, times a constant we will ignore. Given that de Sitter 
space is homogeneous, the diagram is proportional to the infinite spacetime volume 
V,^jd^y^\: 



{z{x,xo) - 1 - ieY 



^Broken by the interactions, however. 
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and we can choose Xq as an arbitrary point. Choosing xq as the north pole = 
(0,/,0,...,0), we get 

/• 4 r— 1 .A2(47r)^^ /• , ,^ cosh^tsin^^ 



(cosh t COS — 1 — ze)^ 
To compute the imaginary part, we can use the well-known formula [9] 

1 1 .TT 



(a; - iO)^ x^ 6 



It follows that 



T A2(4vr),, /■ , cosh^tsin^^ 

ImMo^o = ^ \„: V^ I dtde 



4 



2-41/4 V (coshtcos^-i; 

A^(47r)^^^^ /■ ^a^^g Sin^ (1 + x)3sgn(x, ^) 3) 



2-41/4 y x4 cos 3^ ^(l+a:)2-cos2^ 

where the principal part P regularizes the divergence at x = and sgn(x, 6) is zero 
except in the case (1 +x) cos 6* > 0. It seems clear that the angular integral diverges 
owing to the behavior of the integrand in a neighborhood of 6' = |. 

This is a very strange result indeed, because as we shall see in a moment, this 
is equivalent to a corresponding divergence in the vacuum decay amplitude in the 
tree approximation, which has been reported [10] to be finite in the literature. This 
interesting calculation is reviewed in the Appendix. 

The result we get is 

.A^(47r)^, / s / ivr -ir in \ .A^(47r) 



2-4!/ 



This expression is divergent in the limit when e — )■ O"*", and it has got a finite 
part in this regularization of sorts we are using. Of course we are aware that only 
the limit when e — ?■ O"*" has got any physical sense, so that the divergent parts are 
regularization dependent and just have to be renormalized away. 




Figure 1: Unitarity relation in the vacuum-to- vacuum amplitude to second order. 
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Let us now turn our attention to the computation of the vacuum decay rate. 
We use again the de Sitter homogeneity as we did for the vacuum bubble"^. The 
corresponding integral contains Wightman functions instead of propagators 

|Mo^4p = ^ / d'x^\fldk,hX^)f*^ix,) = 
■ ^ i=i 

= 1 d'x^\W{x,xor = I d'x^\G{x,xor+ 



6\e\ 



(2.5) 

where we are again ignoring the concrete normahzation factors of the two-point 
functions; and we have used that the integral of the propagator over half of the full 
hyperboloid (over t G [0, oo)) is precisely the integral in (2.4), J(e), divided by two. 

The vacuum decay matches exactly the imaginary part of the vacuum energy in 
the regularized theory (as in (1.1)), which is more than is necessary for unitarity, 
which does not necessarily hold for the regularized theory (id est, before taking the 
limit e — 7- 0^ in our case). 

It is worth remarking that the total vacuum decay rate per unit volume and unit 
time interval at tree level appears to be divergent in the e — limit. Given the 
fact that the corresponding tree level amplitude is indeed finite, the divergence is 
entirely due to the integration over the phase space. We are not aware of any such 
divergences, except the ones associated to bremsstrahlung corrections in the external 
legs [11]. It would be interesting to study what happens in the non-conformal case 
as well as to investigate whether this divergences are related to the ones found in 
[12] in a classical context. 

In the Appendix we discuss a different way of computing the vacuum decay rate, 
in a form that results proportional to the covariant four-volume in a way consistent 
with Fermi's Golden Rule; there seems to be an ambiguity as to whether the Rule 
can be applied in arbitrary coordinates. 

Indeed, in [1] some arguments are given for finiteness (after substraction) of the 
massive contribution to the imaginary part of the effective action, and in [3] [5] it is 
similarly argued for the finiteness of the forbidden decay width, which is, as is has 
been already argued for, a closely related quantity through crossing symmetry. 

The result above is quite general, as these unitarity relations are "built-in" in the 
S matrix formalism. In an A^" theory with an arbitrary mass, the relevant identity 
to second order is 

Re j dxdyG{x,yY = j dxdyW{x,yY (2.6) 

^The Wightman functions are not fully de Sitter invariant (not under time reversal). However 
this is enough to repeat the procedure. 



-7- 



Using again the homogeneity of de Sitter, we obtain 



Re j dxG{x,xoT = j dxW{x,xoY (2-7) 

and the relationship between the propagator and the Wightman function G{x,y) = 
T(W{x,y)) allows us to decompose these integrals 

dx G{x, XoT = 1^ + 1- , J dx W{x, Xq)" = /+ + (/")* (2.8) 

where we separate de Sitter space in two regions, future and past of Xq, and those 
are their respective contributions ^. So now our identity looks 

Im/+ = Im/~ (2.9) 

The decomposition of the Feynman propagator is 

G{x, y) = ^{G'^^\x, y) + ia{x, y)D{x, y)), 

where the Hadamard symmetric function G^^^ and the Pauli- Jordan commutator 
function D are real, a{x,y) is the sign of the time ordering of x and y, and D is 
antisymmetric and zero for causally disconnected points, i.e. has the form D{x, y) = 
a{x,y)A{x,y) with A{x,y) symmetric. 



The imaginary part of the nth power of the propagator is then proportional to 
iA(x,y) times asymmetric function (built with odd powers of G^^^ and even powers 
of D), in such a way that the integrand depends only on z{x, xq), and not on the sign 
of the time ordering, so the transformation — )■ —X^ leaves unaltered the value of 
the integral. This proves the unitarity relation. 

There is a formulation of unitarity (the largest time equation)^ which acts of 
position space Feynman diagrams themselves and, as such is most suitable for appli- 
cation in curved spacetimes. This set up is due to Veltman [13], and asserts that for 
any diagram F (xi, . . . , Xn), 

2ReF (xi, . . . ,x„) = - ^ F{xi,...,Xn) 

cuttings 

This formula stems from a representation of the Feynman propagator as 

^Actually the argument x can be in the future, the past or be causally disconnected from 
Xq. We are splitting artificially the spatial region in two pieces (respect to some time parameter) 
and including them in the two causal contributions. This procedure is correct, because both the 
propagator and the Wightman function have the same real value in that region. 

^To be precise, the statement is that in flat space the largest time equation can be shown to 
imply unitarity of the S-matrix. 
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(x, y) = 9{t)A+{x, y) + ^^(-t) A^(x, y) 

such that 

(A+(x,i/))* = A^(x,i/) 

Those conditions are fulfilled in our case. It is always possible to argue that all we are 
doing is to check the largest time equation, without necessarily committing ourselves 
to the thorny [14] issue of unitarity in de Sitter space. 

In reference [3] it is stated that the quantitative importance of the forbidden vac- 
uum decays for massive particles as well as small values of the cosmological constant 
is 

g-'|Am| 

where Am is the violation of energy conservation in the decay process. In fact, 
ml = 1 corresponds roughly to the quintessence scale Hq ~ 10~^^ eV at the present 
period of cosmic evolution. 



3. An imaginary part for the cosmological constant? 

The standard lore on renormalization in Quantum Field Theory [15] includes the 
fact that the cosmological constant is additively renormalized away, the finite part 
being fully undetermined. There is an exception to this, however. In some cases, 
when nontrivial boundary conditions are imposed on the fields (like vanishing elec- 
tromagnetic field in two fixed parallel plates in the classical example) the difference 
between the vacuum energy corresponding to nontrivial boundary conditions and 
the one with "trivial" boundary conditions is computable and in many cases, finite. 
Confer [16] for a clarifying review. This effect is usually known as the Casimir effect 
and has nothing to do with the present work; boundary conditions are kept fixed 
( "trivial" ) in our setting. 

Let us review this in the presence of an external gravitational field. In order to 
do that, we shall consider a simple model of a scalar quantum field 



(the terms quadratic in the curvature are necessary for renormalization). We do not 
believe that our physical conclusions depend upon the details of the model. Standard 
computation of the one-loop effective potential [17] yields the result 

VeS! = Ao + \iR¥ - azR^ - a^R^'^'^R^,,^ - a2R^'' R^u + \ml^l + ^0^+ 
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Defining renormalized quantities in the MS scheme 

\ n — A 
(where A is dimensionless) leads to 

n .A h 

/3a = -4A + — — 

The physical Cosmological Constant (CC) ^ is given by 

^phys = A'^A 

so that 

R - 4 

P^phys — 327]-2"^ 

Let us define the cosmological constaw^ in the MS scheme as 

A(A)^K/K0 = O) 



Assuming 



<l>=0 



0, that is, absence of spontaneous symmetry breaking, 



the CC obeys the renormalization group (RG) equation 

A successfully employed technique to compute the effective potential in many cases 
[8] stems precisely from using this equation in a recursive way. Feeding the RG 
equation with the MS one loop result 



T m H -m log —— 

327r2 V2 4; 647r2 ^ 



^This is entirely analogous to 't Hooft's [18] beta function for the dimensionless mass 

Xm 



327r2 

in such a way that the physical mass 

rriphys = rnfi 

and 

n Xnipfiyg 

P^^pHys = 327^2 

^Of course the CC is put by hand equal to zero if we were to use Coleman- Weinberg's [19] 
renormalization condition 

^e//(0 = O)=O 
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and using the well-known one-loop beta functions and anomalous dimensions ^ yields 
the dependence of the two loop piece: 



A2 



256^ V2 




h 



Xra^ log h 



hXm^ 



) 



log jji— 



10247r4 



h 



Am^ {logiif + f{gi) 



5127r4 



where gi stands for all coupling constants not including the CC itself, but including 
the this piece is not determined by the RG equations. 

No imaginary part appears ever by this procedure. Some ambiguity remains in 
the finite part independent of the coupling constants, which is usually fixed by the 
renormalization conditions, that is, renormalized to zero ([8]). 

The aim of the present work is to point out, following the logic of [1] and the pre- 
vious section that generically owing to lacking of translation invariance, the matter 
vacuum is unstable (that is, it can decay to physical particles with a certain com- 
putable width) so that by unitarity the vacuum energy must have a unambiguous 
finite imaginary part. Assuming analyticity in the coupling constants, this in turn 
puts restrictions on the possible real part of the cosmological constant through the 
dispersion relations as a consequence of Cauchy's theorem. This explicitly contradicts 
previous claims in the literature [20]. 

Another consequence of all this is that the usual analysis of quantum field theory 
in curved space-times is incomplete without a careful self-consistent consideration of 
the back-reaction problem, owing precisely to this imaginary part. We shall comment 
on this point in the conclusions. 



^To wit: 




180 167r2 



(3.1) 



4. A model with unstable vacuum in flat space 



We have just witnessed a few paragraphs ago that exphcit computations even in the 
simplest of background spacetimes, such as constant curvature de Sitter or anti de 
Sitter spaces, are quite involved [3], and besides, the physical interpretation of the 
divergences found is not clear. 

This is the rationale for first performing the complete analysis in a much simpler 
context in flat space, a toy model of sorts. It is clear that the main new ingredient of 
a curved space in this context is the non-conservation of four momentum. This can 
be achieved through an explicit time-dependent interaction term in fiat space which 
violates energy-conservation. In the appendix the mapping between a subclass of 
scalar models in Friedmann spacetimes and fiat models with spacetime dependent 
coupling constants is spelled out. 

Let be the lagrangian 

^ = (4.1) 

with a non- homogeneous coupling A(t) = |(l + cos?7t). With this choice, energy 
is no longer conserved and a fixed amount {rj) can be created or destroyed at any 
vertex. In the limit 77 —t- we recover an standard X(j)^ theory. 

With this interaction, the second order vacuum-to-vacuum amplitude reads: 



^ = Texp j A(x°) : (f){xf : c/x| (4.2) 



(OI^^'^IO) = y \{x')X{y')G{x, yf dxdy = 

_ .(27r)W f,o,of dp. \ , o^n , oJHpt) e^^ri^'-y"^ 



J l^^^l (1 + cosr/x°)(l + cosr/?/ 



Dpi Dp^ Dp^ 



(4.3) 



where we abbreviate by Dp the p"^ — m? + ie denominators, {dpi} and pt mean dpi 
and ^jPi; V is an infinite volume factor. 

The time dependent coupling can be expressed as proportional to 1 + ^e^^^° + 
ig-irjx"^ so the "energy conservation" factors are now different 

(0|5(2)|o) = 



8-3! J \{2nr ! Dp,Dp,Dp, {'^P^^ ^ 4 
.(27r)4yTA2 ^ ^ 1^ ^2 



i ^-^^^ (T234(0, m,m,m) + -T234(?7 ,m,m,m) ] (4.4) 
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so the creation of energy at one vertex should be compensated in the other. T is 
an infinite time factor. The standard integral T234 is well known [21] in general 
dimension: 



7234(p^m2,m3,m4) 



d"A;2rf"A;3d"A;4 S'^ih + ^2 + ^3 - p) 



(27r)3'^ {k1 -ml + ie) {kl - m| + ie) [kj -ml + ie) 

(4.5) 

and corresponds to the self energy "setting sun" diagram S(?7^) in a theory. 

Notice that we have assumed that 77 7^ 0, so the regular theory result is not 
recovered when 77 —j- in the last formula, i.e. the limit is discontinuous. 

The corresponding vacuum decay rate can be calculated as the square of the 
T-matrix up to first order, T^^^ ~ I dxX{x'^) : 4){x)^ : 

ro^3 = (o|r(^)r«t|o) = 

= dxdy{0\ : 0(a;)3 :: 0(y)3 ■ \0) = \{x'>)X{y')W{x,yy' dxdy (4.6) 
We ca use now the expression for the Wightman function of the scalar field 

W{x,y) = I ^(27r)5+(p^-m^)e-^^(^-^) = J ^^/''^''-'^e-^'''^^'-''^ (4-7) 
and then 

(27r) VTA^ f (27r)5+(p2 - m'^)dp. 



3!4 J (27r) 
(27r) VTA^ f (27r)5+(j92 - m^)dp, 



3!16 



{2TTyVTy f dpidp2 
3116 



where -£"12 = a/ (Pi + P2)^ + 'n^-'^- From the three delta function factors for te energy, 
only the one with positive increment of energy contributes, for kinematical reasons. 

This expression is just the standard three-body phase space factor 
_ {2nrVTX' [ dE.dE, 

- 3!16 j 4(2^0^^^^^ -r?)%-i?T ) (4-9) 

where iJmhf = -|_ + ^ [p^ ± -\- m? correspond to parallel and anti-parallel 
configurations of momenta. There is also an implicit kinematic factor 6{r] — 3m) 
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indicating the minimun energy necessary to create 3 particles. 



The final expression for this rate is 

2,2 

^0^3 = - — / — ^HE,r]) 4.10 

with 

^y (E^ — m?){rf + m? — 2Eri){rf — 3m^ — 2riE) 
nE,ri) = ^—^ — (4.11) 

This rate is proportional to the decay rate of a particle with mass t], TM=ri- We can 
give he approximate result of this integral if rj is very close to or much bigger than 
3m: 

2 „ 2 
77 —3m 

/ dEI{E,r]) ^2V3m^e'^ T] = 3m{l + e) (4.12) 

J m 

2 2 

—3m 2 

(iEJ(E, 77) ~ ^ , if r/ > m (4.13) 



From the application of Cutkosky's rules, the "setting sun" self-energy diagram 
has an imaginary part given by the corresponding decay rate, ImS(?7^) = TM=ri^ 
which is based in the identity [21] 

2 2 

77 —3m 

/ dE/(E,r/) = -^Im((27r)3«T234)'" (4.14) 

Jm 47r5 „=4 

Since the processes in our model are directly related to the aforementioned dia- 
grams, we can establish easily the unitarity relation 

-2Re(0|^(2)|0) = ro^3 (4.15) 

In the previous calculation, we avoided to deal with the tadpole diagrams by 
taking the interaction operator to be normal ordered. If we choose to include those 
extra diagrams, the unitarity relation still holds. The tadpole contribution in (4.15) 
is then: 

-2Re(0|5(2)|0) = ro^3 + Tq^i (4.16) 
Here the left member acquires a new term which is 

(0|521|0) = -\ f dxdy\{x')\{y'^)G{x,x)G{y,y)G{x,y) (4.17) 
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The closed loops G{x, x) and G{y, y) are constants, and they are always real (at least 
in dimensional regularization) and so, not very relevant for the calculation. 

i 2 /" dxdydp on/o o q e*^*^^"^) 



(0|4^>) = -^Gg / ^g^(2 + 2cosr^x°)(2 + 2cosV 



Dp 



^ ' ^Gl / dxVdp°(2 + 2cosr7,T°)(2 + 2cosr7y°^ 



8-16 V • ™ ^ (p0)2 _ ^2 + 

'''^'''''Ol(±.^^^\ (4.18, 



8-16 ° 7]'^ — m? + ie^ 

In the other hand, the total decay rate ro-s>i is 

ro^i = ^ E / dxdyX{x')\{y'){m^f\^){myf\^) = [ dxdyX{x')X{y')W{x,y) 



3! 

1— particle 



l^Gl I ^^^^e^P('^--)(2 + 2cosr^x°)(2 + 2cosr^y°)e 



0^ iE^ix"-y^) 



64 

= ^^^^oi ^^^^ - + + ^^-^9) 
We have assumed W{x, x) = G{x, x). 

These new contributions to (4.15) should match each other. It is easy to check 
that they actually do by applying the Weierstrass theorem 

^ ^ {6{r] - m) + 6{r] + m)) + . . . (4.20) 



7]^ — m? + ie 2m 

where the dots indicate imaginary contributions. No infinities nor ambiguities have 
been encountered in our toy model. 
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5. Kinetic equations 



The effects shown in the previous section, allowed because of the non-conservation 
of energy, led to the production of particles in the initially empty space. In our very 
simplified model there are no conserved quantum numbers; but in general the set of 
produced particles must enjoy the quantum numbers (id est, charges) of the vacuum. 
An important physical question is towards which final state this instability leads to? 
Again, this question has many faces. In the general case in which a gravitational ex- 
ternal field is present, the backreaction is surely important, but difficult to compute. 
It is not even clear that the usual procedure of solving again Einstein's equations 
with a second member given by the expectation value of some energy- momentum 
tensor [22] would be good enough for our purposes. 

A second facet of this problem is the evolution of the particle (in general, con- 
served quantities) density. We can study this phenomenon by considering a spacetime 
box V X T, with an initial number of particles (in our simplified model there is only 



one type of particles) A^, do that the initial density is n _ ^ 



N 10 



It is well-known [23] that the transition rate dT {a — )■ /3) = '-^^ of a given 
process in which Na particles in the initial state evolve into Nf^ particles in the final 
state depends on the three-space volume, V, as V^'^". This clearly means that the 
vacuum decay terms in Boltzmann's equation (or whatever improvement thereof) 
will clearly dominate, because they are the only ones which are extensive, id est, 
proportional to the ordinary volume. Let us build up a simplified model for this 
equation. The essential thing is to capture the volume as well as the power of the 
distribution function itself. Taking into account the vacuum decay and absorption 
only, the balance equation in our fiducial spacetime box is of the schematic form 

dt ~ 

Here A^p dp represents the total number of particles with momentum between p 
and p + (ip in the fiducial volume V, and T is our fiducial time. The constructs 
and F^q denote the amplitude for vacuum decay or annihilation to the vacuum 
when precisely one of the particles has momentum p. The transition rate is simply 
related to the transition probability dP{a (5) by dV{a (3) = General 
arguments implicate 

dT{a /3) = {27rf''--^V^-^-\Mp^\^6\p^-pfs)d/3 

where the matrix M is intimately related to the matrix S: 



^"^A finer description would be provided by the Wigner function, in wliicli both momentum 
and position distributions are correlated, to the extent that this is compatible with Heisenberg's 
principle. This construct obeys suitable generalizations of Boltzmann's equation. 
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and 

c?/3 = J]^ dpi 



|2 



With this proviso, the mass dimension of \M\ is 8 — 3(iVa + N/s)- 

It is now clear that the increase in the density of particles of momentum p out 
of the vacuum in the fiducial volume V during the time interval T is 

rP3 = J dro3ip,p'^,p',)dp'^dp',ii + Np)ii + Np,)ii + Np,). 

The decrease in the density of particles of momentum p due to annihilation to the 
vacuum (id est, the reverse process) in the same fiducial volume and time interval is 
proportional to the existing density cube of particles: 



rP = 

^ 30 — 



= J NpNp^Np^ dp2dps dr-io (p,P2,P3) 
The total (integrated) widths are given by 



2 rj —3m 2 

where the last equality holds approximately when rj — 3m(l + e). 

Tao = J dpidp2dp3Np^Np^Np^dT3o{pi,p2,P3) = 

^^(27r)^ f fipi 0?p2 rfps 3 , , ^ f.o\ 

^^GV^J ^-^-^-21:2^21^^ (P1 + P2 + P3) (5.2) 

A toy model that embodies some of the characteristics of our integro-differential 
equation is 



dt V M\ 

where C is a constant and we have saturated all dimensions with an average mass 
scale, M. It is clear that in any such model the density of particles will rapidly grow 
as 

n ~ CMH 

until it becomes so big that 
which happens in a characteristic time 

1 

CM 
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6. Conclusions 



In this paper the issue of vacuum instabihty in the presence of a cosmological constant 
has been discussed. Some computations for conformally invariant matter in de Sitter 
space were done, and an ambiguity with Fermi's Golden Rule discussed in some 
detail. 

A flat space model with time-dependent coupling constant, that captures some of 
the aspects of the physically interesting situation has been introduced, and analyzed 
to one-loop order. No ambiguities have been encountered in this piece of work. 

It is somewhat hazardous to transport all lessons from this model to the quantum 
fields interacting in an external gravitational background, but it is clear at any rate 
that vacuum instability is a generic phenomenon for any quantum field in the presence 
of such an external gravitational background. Any attempt to forbid it by using 
clever quantum numbers seems doomed by the fact that it would always be possible 
to create particle-antiparticle pairs with vacuum quantum numbers. Now there are 
not so many logical possibilities. 

One that immediately comes to mind is that unitarity can be violated. We have 
pointed out already in section 2 that what is essential here is not the full fledged 
notion of unitarity, but only the milder consequence of the definition of Feynman's 
propagator that is usually known as the largest time equation. In flat spacetime 
this is known to implicate unitarity, but perhaps in a general background this need 
not be so; nevertheless the largest time equation itself seems unavoidable as long as 
Feynman's propagators are used. It could be that ordinary perturbation theory gets 
modified in this context; careful computations [3] seem to indicate that this is not 
the case, although the final word has not been said . 

If unitarity in this sense holds true, then, as we have argued, the cosmological 
constant gets a nonvanishing imaginary part, which points to the existence of some 
instability. Now what could be the endpoint of such an instability? It is clear that 
the only spacetime in which this instability is absent is flat space-time. 

It seems difficult at any rate to resolve this issue in the hybrid framework [24] 
in which the quantum expectation value of the energy-momentum tensor is used as 
a source for Einstein's equations. It is likely that the gravitational field has to be 
treated in equal terms with all other fields in order for a self-consistent solution to 
exist. 

Detailed calculations in this framework are necessary but difficult. We hope to 
report on them in the near future. 
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A. Alternative computation of the vacuum decay rate. 

Let us review the calculation in [10]. In the conformally flat patch 

ds^ = — (du^ — dx^) 
u 

for a conformally coupled scalar fleld we can expand the free fleld expansion in 
term of the modes (we represent by A; = |k|). 

which are equivalent to choosing the Euclidean or Bunch-Davies vacuum^^. The 
normalization is such that [01,03] = (27r)^(5i2. This expansion is the same for the 
interacting fleld in the interaction picture, so the amplitude — > kik2k3k4 is simply 

dndx(kik2k3k4|-0(n,x)|O) = A(27r)^2 / dx JJ /^^ (n, x) = 

^ i=i 

= X [ /"^^ e''"g-^'e-g--'-' = A(27r)3 [ , e-^^^- ^'^(V k,) 

(A.l) 

The total decay rate 

We consider that the inflnite quantity 

K = (27r)35(0) 

represents the volume of 3-space. The following changes of variables are considered: 

u = I e i , 2T = ti + t2 

v = le~'^, T = ti-t2] (A. 2) 



^^A conformally coupled scalar field is massless and the value of the curvature coupling is 
^ = ^i^I^' 'v^^here n is the dimension of spacetime. In de Sitter space, where the curvature 
is constant, R = — "^"a'"'"^ , this is equivalent to a minimally coupled = 0) scalar field with 
^iri^^ = '-^^^jiT^. This mass is in the complementary series, with a ifi — pa- 
ramctcr (m^P = fi^ + ^"^ ' ) , and its corresponding two-point functions arc particularly simple 

^^For the n-dimcnsional case, the appropriate modes are /k(u,x) = ^^*-'^"^^^^ e'*^"e''^^ 
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yielding 



4!(27r)n6 



/" dT dr H ^6iJ2 k^e-'T exp 1 2tl ^ h e"^^ sinh ^ 1 

i=i i=i y i u 



Now the integral over the momenta includes a term which looks singular 

1 



/ 



dki(£k;2(ik3 



A;iA;2A;3|ki + k2 + ka 



The point is that |ki + k2 +k3| can well vanish even if all /cj 7^ 0. Expanding around 
the seemingly singular region, it appears that the singularity is integrable. 

At any rate, the integral diverges for high momenta, unless the exponential 
decreases, which happens only when 

Im ^ fcj > 

i 

Nevertheless, it can be defined as a Gamma function by analytic continuation, 
so that let us carry on. 

In the reference we are annotating unity is introduced in the form 

poo 4 

•^0 i=i 

and the momenta are normalized by kj = Ky^, so that using 

i=l i i 



3 
4" 



yields 



P = ^ / dTdrdKK'e-'-t exp | ^e"? sinh 
= j dTe^T Vc (^ 4g^4^^g^-)6 ^'^'^^ J exp(i«; sinh rj) 

T 

where k, — 2Kle~T and 11—21 particle production rate per unit volume is 



— dKK'J d77exp(msinhr^) = ^^^j^ dnnHK^iK) 



This formula is very appealing physically; it gives a finite result for a tree-level cross 
section, which is a good thing, because the only known place in which divergent cross 
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sections appear at tree level is in bremsstrahlung effect (correction to external legs 
by emission of massless particles), and this is not our case. 

This means that whereas the calculation in section 2 was proportional to the 
covariant spacetime volume, V4 (times a divergent expression), Higuchi's is propor- 
tional to Vc times the integral over dTe^^ times a finite expression. If we identify 

1/4 ~ 1/3 y dTe^T 

both calculations are inconsistent. But this identification is not compulsory. As 
has been already pointed out in [3] it is not clear in which coordinates the particle 
production per unit volume per unit time is to be defined. 

Ay any rate, by consistency, in case the computation of [10] is preferred, the 
imaginary part of the free energy should also be finite. 



B. Relationship between time-dependent coupling constant 
and cosmological backgrounds 

Let us look for the cases in which a Minkowski space lagrangian with time-dependent 
coupling constants such as 

is equivalent to a scalar field in a curved, conformally flat gravitational background 

ds"^ = aitY'ri^ydx^dx" 

namely 



^rd,4,dA ' '^4? - li) = y (*' - (*')') - o'^*' - a^l'P" 



This is always possible, up to a total derivative, with the identifications 

ip = a(l) 
\{t) = \a{t) 

m\t) = m'a^ - - (B.l) 
a 

It would also be interesting to find a mapping for the different vacuum states as 



well. 



^^In the last formula the integral representation for the Bessel function Ki^{z) has been used. 
This is for -7r/2 < Arg z < 7r/2 

00 poo |'00 + ?.7r/2 



/ i7r/2 

which seems to require that Imrj ~ 7r/2, which is not the case, but again the integral can be 
interpreted as the analytical continuation of a convergent one. 
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